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The NLO QCD calculation for the photoproduction of the isolated photon with a large pt at the HERA ep collider is presented. 
The single resolved photon contribution and the QCD corrections of order a s to the Born term are consistently included. The 
NNLO contributions, the box and the double resolved photon subprocesses, are sizeable and are taken into account in addition. The 
importance of the isolation cut, as well as the influence of other experimental cuts on the pj- and rj-y (the final photon rapidity) 
distributions are discussed in detail. The investigation of the renormalization scale dependence is performed in order to estimate the 
size of missing higher order QCD corrections. Results are compared with experimental data and with the prediction of a different 
NLO calculation. 



1 Introduction 

The production of the prompt photon with large trans- 
verse momentum, p?, in the ep collision is considered. 
Such reaction is dominated by events with almost real 
photons mediating the ep interaction, Q 2 sa 0, so in prac- 
tice we deal with the photoproduction of the prompt pho- 
ton. The other name for such process is Deep Inelastic 
Compton (DIC) scattering (although Q 2 as 0, the scat- 
tering is "deep inelastic" due to a large transverse mo- 
mentum of the final photon) . The photon emitted by the 
electron may interact with the proton partons directly or 
as a resolved one. Analogously, the observed final pho- 
ton may arise directly from hard partonic subprocesses or 
from fragmentation processes, where a quark or a gluon 
decays into the photon. 

The importance of the DIC process in the ep colli- 
sion for testing the Parton Model and then the Quantum 
Chromodynamics was studied previously by many au- 
thors a a. Measurements were performed at the HERA 
ep collider by the ZEUS groupB-feJ, also the HI Col- 
laboration has presented preliminary results 113. In these 
experiments only events with isolated photons were in- 
cluded in the analysis, i.e. with a restriction imposed on 
the hadronic energy detected close to the photon. The 
corresponding cross sections for the photoproduction of 
an isolated photon and of an isolated photon plus jet_werc 
calculated in QCD in next-to-leading order (NLO)EJt£|. 
There also exists analogous calculation for the large-Q 2 
ep collision (DIS events) i3. 

In this paper the results of the NLO QCD calculation 
for the DIC process with an isolated photon at the HERA 
ep collider are presented. We consider the parton distri- 
butions in the photon and parton fragmentation into the 
photon as quantities of o-cder_a em . Our approach differs 
from the NLO approach EaTla by set of subprocesses in- 



cluded in the analysis. The comparison of our predictions 
with the NLO results obtained by L.E. Gordon (LG) Ej 
is presented for cross sections with kinerpatical cuts as in 
the ZEUS Collaboration measurements L3. 

The present analysis is the final, .much extended and 
improved version of the previous onet-3. We show results 
for non-isolated final photon, and we study the influence 
of the isolation cut on the production rate of the photon. 
The role of other specific cuts applied by the ZEUS CoL 
laboration are discussed and the comparison with dataE3 
is made. We emphasize the importance of the box di- 
agram 7<? — > 7(7, being the higher order process, in de- 
scription of the data. 

We study the renormalization scale dependence of 
the cross section in order to estimate the size of missing 
higher order (NNLO or higher) QCD corrections. The 
NLO results for the photoproduction of the isolated 7 
are compared to the leading logarithm (LL) ones, and in 
addition the LL predictions for isolated j+jet final state 
are presented. 

In the recent ZEUS analysis of the prompt photon 
plus jet productionEEl the intrinsic transverse momentum 
of partons in the proton was included in Monte Carlo 
simulations to improve agreement between data and pre- 
dictions. This momentum is not included in our calcula- 
tions. 

We start with discussion on the choice of relevant 
diagrams defining our NLO approach to the DIC process 
(sec. ||). The isolation of the photon is described in sec. |[ 
and the equivalent photon approximation in sec. [|. In 
sees. H and ^ the results of numerical calculations are 
presented and compared with data E3 and other NLO 
predictions Il3. In sec. [?] we show LL predictions for the 
photon plus jet production. Finally, sec. || summarizes 
our results. 



Figure 1: The Born contribution. 

2 The NLO calculation for •yp — > 7X Deep In- 
elastic Compton scattering 

2.1 General discussion 

We start by describing processes which are (should be?) 
included in the NLO QCD calculations of the cross sec- 
tion for the DIC process, 

IP -> lX, (1) 

where the final photon is produced with large transverse 
momentum, pt ^> Aqczj- Although we will consider the 
process (|l]), the problem which we touch upon is more 
general - it is related to different approaches to NLO 
calculations of cross sections lox. hadronic processes in- 
volving resolved photons, see&LD and for more detailed 
discussion E9. 

The Born level contribution to the cross section for 
process (|l|), i.e. the lowest order in the strong coupling 
a s term, arises from the Compton process on the quark 

(%• 0): 

11 -> 7<7- (2) 

It gives the \a 2 m ] order contributions to the partonic 
cross section q At the same a 2 m order it contributes 
to the hadronic cross section for the process "fp — > 7A. 

The Parton Model (PM) prediction for the DIC pro- 
cess (|l|), which applies for xt — 2par/VS ~ C(l), relies 
solely on the Born contribution (0)u, namely: 

da iv->ix =J2 I dx p q p {x p )da^ q ^\ (3) 

where q p is the quark density in the proton and d<7 lv ^ lX 
(da' lq ^ lq ) stands for the hadronic (partonic) cross sec- 
tion. In the QCD improved PM the cross section is 
given by @, however with scale dependent quark den- 
sities. For semihard processes, where it < 1, the pre- 
diction based on the process (g) only is not a sufficient 
approximation, and one should also consider the contri- 
butions corresponding to the collinear showers, involving 
hadronic- like interactions of the photon(s). There are 
two classes of such contributions: single resolved with 
resolved initial or final photon, and double resolved with 

a We denote the order of partonic subproccsses using square 
brackets. 



Figure 2: Examples of single resolved processes: a) the resolved 
initial photon, b) the resolved final photon. 




Figure 3: An example of a double resolved photon process. 

both the initial and the final photon resolved (figs. |[ 
||) . They correspond to partonic cross sections of orders 
[a em a s ] (single resolved) and [a 2 s ] (double resolved). If 
one takes into account that partonic densities in the pho- 
ton and the parton fragmentation into the photon are 
of order ~ Qt em , then the contributions to the hadronic 
cross section from these resolved photon processes are 
a 2 m a s and a 2 m a 2 , respectively. Both single and double 
resolved contributions are includcd-,in the standard LL 
QCD analyses of the DIC process □!§□. 

To obtain the NLO QCD predictions for the pro- 
cess (fy) the a s corrections to the lowest order pro- 
cess (0) hasz^jkj-be calculated leading to terms of or- 
der oq m a s □•BEilS (fig. ||). In these a 2 em a s contribu- 
tions there are collinear singularities to be subtracted 
and shifted into corresponding quark densities or frag- 
mentation functions. This way the single resolved photon 
contribution appears in the calculation of the a s correc- 
tions to the Born process. It is worth noticing that in 
the NLO expression for the cross section there are no 
collinear singularities which would lead to the double re- 
solved photon contributions. It indicates that taking into 
account [a 2 ] subprocesses, associated with both the ini- 
tial and final photons resolved, goes beyond the accuracy 
of the NLO calculation. This will be consistent within the 




Figure 4: Examples of the virtual gluon and real gluon a s correc- 
tions to the Born contribution. 
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NNLO approach, where a 2 correction to the Born term 
and a s correction to the single resolved terms should be 
included, all giving the same oQ m OQ order contribution 
to the hadronic cross sections. 

Thfi other set of diagrams is considered by some au- 
thors Ed in the NLO approach to DIC process (|l|), 
due to their different way of counting the order of par- 
ton densities in the photon (and the parton fragmen- 
tation into the photon). This approach, which we will 
call "l/a s " approach, is motivated by large logarithms 
of Q 2 in the existing already in the PM. By ex- 
pressing ln(Q 2 / Aq CD ) as ~ one treats the par- 
ton -densities in photon as proportional to a em /a s (see 
e.g.y^ffitrlla). By applying this method to the DIC 
process, we see that the single resolved photon contribu- 
tion to the hadronic cross section for 775 — > ^X becomes 
of the same order as the Born term, namely 



[a em a s ] <g) 1 = a 2 . 



(4) 



The same is also observed for the double resolved photon 
contribution 



® at 



(5) 



We see that by such counting, the same a 2 m order 
contributions to the hadronic cross section are given by 
the direct Born process, single and double resolved pho- 
ton processes although they correspond to quite different 
final states (observe a lack of the remnant of the pho- 
ton in the direct process). Moreover, they constitute the 
lowest order (in the strong coupling constant) term in 
the perturbative expansion, actually the zeroth order, so 
the direct dependence of the cross section on the strong 
coupling constant is absent. If one takes into account 
that some of these terms correspond to the hard pro- 
cesses involving gluons, the lack of terms proportional to 
a s coupling in the cross section seems to be contrary to 
intuition. 

In the "l/a s " approach beside the a s correction to 
the Born cross section, the a s corrections to the single 
and to the double resolved photon contributions are in- 
cluded in the NLO calculation._siiLce all of them give 



terms of the same order, a. 2 



To summarize, the first approach starts with one ba- 
sic, direct subprocess as in the PM (eq. |), while the sec- 
ond one with three different types of subprocesses (as in 
the standard LL calculation). Obviously, some of NNLO 
terms in the first method belong to the NLO terms in 
the second one. 

In this paper we apply the first type of NLO approach 
to the DIC process, however with some important NNLO 
terms additionally included. A comparison between-pps 
results and the res ults based on the other approach 0Jil3 
is discussed in sec. 3.3. 



Figure 5: The box diagram. 



2.2 The cross section 

Below we describe our approach to the DIC process, 
where the parton densities in the photon and the par- 
ton fragmentation into the photon are treated as ~ a em . 

In the NLO QCD calculation of the DIC process we 
take into account the following subprocesses: 

• the Born contribution ^ (fig. |l|); 

• the finite a s corrections to the Born diagram (so called 
K-term) from virtual gluon exchange, real gluon emission 
(fig. ||), and the process 73 — > qqj; 

• two types of single resolved photon contributions, with 
resolved initial and final photons (fig. |^). 

Besides the above full NLO set, we will include two 
terms of order a 2 m a 2 (formally from the NNLO set): 
the double resolved contributions (fig. |j) and the direct 
diagram (boxl 717 — > jgc3 (fig. ^|), since they were found 
to be large Uu. 

The cross section for the 7p — > jX scattering has the 
following form: 



2n 2 s 



-K h 



/ ~J j dX 1 j dx fa/~f( x liQ 2 )fb/p(x,Q 2 ) 

abc 



■D l/c (z,Q 2 )E 7 



d 3 a a 



b— >cd 



(6) 



The first term is the K-term describing the finite a s cor- 
rections to the Born process, and the second one stands 
for the sum over all other contributions (including the 
Born contribution). The / a / 7 (fb/p) is a a (6)-parton 
distribution in the photon (proton) while the D 1 / c is 
a c-parton fragmentation function. For the direct ini- 
tial (final) photon, where a — 7 (c = 7), we take 
/ / 7 = 8{x 1 — 1) (D y / C — 5(z — 1)) (the Born contri- 
bution is obtained for a = 7, b = q and = 7). The 
variables and z stand for the fraction of the initial 

photon, proton, and c-parton momenta taken by the a- 
parton, 6-parton, and the final photon, respectively. The 
renormalization scale is assumed equal to the factoriza- 
tion scale and is denoted as Q. 
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3 The isolation 

In order to observe photons originating from a hard sub- 
process one should reduce backgrounds, mainly from i s 
and 7's radiated from final state hadrons. To achieve 
this, isolation cuts on the observed photon arc introduced 
in experimental analyses. The isolation cuts are defined 
by demanding that the sum of hadronic energy within 
a cone of radius R around the final photon, where the 
radius R is defined in the rapidity and azimuthal angle 
space (see eq. [l^ in Appendix), should be smaller than 
the final photon energy multiplied by a small parameter 
e: 



hadrons 



(7) 



The simplest way to calculate the differential cross 
section for an isolated photon, dai so i, is to calculate 
the difference of a non-isolated differential_crpss_section, 
da n0 n-isoi, and a subtraction term, da su b CJc^lij: 



duisoi = da n . 



on—isol 



da 



sub • 



(8) 



The subtraction term corresponds to cuts opposite to the 
isolation cuts, i.e. hadrons with the total energy higher 
than the photon energy multiplied by e should appear 
within a cone of radius R around the final photon. 

The isolation cuts are imposed only when calculat- 
ing the K-term, and in contributions involving fragmen- 
tation function (resolved final photon). Other contribu- 
tions arise from 2 — > 2 subprocesses with direct final 
photon that is isolated by definition. 

In the analysis we apply the subtraction method with 
the-subtraction term calculated in an approximate way, 
seeEj'E-3 for details. The approximation bases on the as- 
sumption that an angle S between the final photon and 
a parton inside the cone of radius R is small. It al- 
lows to simplify considerably the calculations and leads 
to compact analytical expressions for all relevant matrix 
elements involved in da su b- Note that in this approxi- 
mation the angle S is simply proportional to the radius 
R: 5 = R/cosh(rj^). The above small-<5 approximation is 
used only in calculation of the i^-term in the subtraction 
cross section da su b, for the results see Appendix. Other 
contributions to da su i,, as well as da non -i SO i and all LL 
expressions, are obtained in an exact way. 

It is worth mentioning that there is an ongoing dis- 
cussion whether the conventional factorization breaks 
down, and whether the cross section is an infrared safe 
quantity for isolated photon photoproduction_iiL n e + e _ 
collisions (also for hadron-hadron reactions) EtE3. In 
principle these questions could as well occur for the pho- 
toproduction of isolated photons in ep collisions. How- 
ever we do not deal with this problem because it arises 



from 2-^3 subprocesses in which a final quark fragmen- 
tates into a photon. We checked this explicitly and found 
that all singularities in da su b are canceled or factorized, 
as in da non ^i SO i. Therefore the considered . by us cross 
section doi SO i is well defined (see also oE£ni3) . 



4 The equivalent photon approximation 

We consider the production of photons with large trans- 
verse momentum, px Aqcd, in the ep scattering, 
ep — * e-fX, at the HERA collider. This reaction is dom- 
inated by photoproduction events, i.e. the electron is 
scattered at a small angle and the mediating photon is 
almost real, Q 2 » 0. The cross section for such processes 
can be calculated using the equivalent photon (Williams- 
Weizsacker) approximationca which relates the differen- 
tial cross section for ep collision to the differential cross 
section for jp collision. For the DIC scattering the ap- 
proximation has the following form: 



da 



ep — >e~yX 



= / G l/e (y)da^ x dy, 



(9) 



where y is (in the laboratory frame) a fraction of the 
initial electron energy taken by the photon. 

We apply the equivalent photon approximation and 
take ±he (real) photon distribution in the electron in the 
formEj: 



2tt 



)2 

■ max 



(i-v). 



m^y z 

? ? 
m e y 

~ Q 2 

max 



)}, (io) 



with m e being the electron mass. In the numerical cal- 
culations we assume Q 2 max as 1 GeV 2 what is a typical 
value for the recent photoproduction measurements at 
the HERA collider. 



5 The results and comparison with data 

The results for the non-isolated and isolated photon cross 
sections are obtained in NLO a ccur acy with additional 
NNLO terms, as discussed in sec. b.2l We take the J3ERA 



collider energies: E e =27.5 GeV and E p =820 GeVB, and 
we consider the pt range of the final photon between 5 
and 20 GeV (xt from 0.03 to 0.13). The calculations are 
performed in MS scheme with a hard (renormalization, 
factorization) scale Q equal pr- Also Q = pt/2 and 2px 
are used to study the dependence of the results on the 
choice of Q. We neglect the quark masses and assume the 
number of active flavors to be Nf=A (and for comparison 
also Nf=3 and 5). The two- loop coupling constant a s is 
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used in the form 

( m = 4lT n - 2(31 HHQV^qcd)] , 

asW > PoHQV^ QCD Y PI ln(0 2 /A| ra ) J 

(11) 

(Pa = 11 - 2/3AT/ and ft = 51 - 19/3JV/), with 
A QC£) =0.365, 0.320 and 0.220 GeV for N f =3, 4 and 5, 
respectively, as fitted by us to the experimental value of 
a s {M z ) = 0.1177B. The tJg CD = 0.120 GeV for N f =4 
was taken in one-loop a s when calculating the cross sec- 
tion in LL accuracy. 

We use the GRV parametrizalions of the proton 
structure function (NLO andJjO)!^, the photon struc- 
ture function (NLD and LO) E3, and the fragmentation 
function (NLO) L3. FoLComparisan other parametrizfL: 
tions arejilso used: DO0, ACFGpI, CTEqH MRSTB 
and GS0. 

As the reference we take the GRV NLO set of parton 
distributions Era, N f = 4, A QCD = 320 GeV and Q = 
Pt- 

5.1 Non-isolated versus isolated photon cross section 

The pt distribution for the produced final photon with- 
out any cut is presented in fig. ^ where the NLO results 
and separately the Born term (with NLO parton den- 
sities) are shown. The cross section decreases by three 
orders of magnitude when pt increases from 4 GeV to 
20 GeV, and obviously the most important contribution 
is coming from the lowest p? region. The subprocesses 
other than the Born one give all together contribution al- 
most two times larger than the cross section for the Born 
subprocess alone. 

The importance of particular contributions to the 
non-isolated cross section integrated over 5 GeV< pr < 
10 GeV is illustrated in Table 1 (the first line). The 
total NLO cross section is equal to 226 pb, with indi- 
vidual contributions equal to: Born = 36.3%, single 
resolved = 35.1%, double resolved = 18.7%, box = 
6.2%, K-term=3.9%. We see that the single resolved 
photon processes give a contribution comparable to the 
Born term. Also the double resolved photon processes are 
important. It is worth noticing that the overall double 
resolved photon cross section is build from many, rela- 
tively small, individual terms. The direct box diagram 
(7.9 75) gives 17% of the Born (75 — > jq) contribu- 
tion. Relatively large box contribution (although being 
[o^]) is such partially due to large gluonic content of the 
proton at small x„. 

Next, in fig. p] we compare the differential cross sec- 
tion da/drjy for the non-isolated photon with correspond- 
ing predictions for the isolated photon using various val- 
ues of the isolated cone variables (e, R) []. The isola- 

b The positive rapidity is pointed in the proton direction. 



tion cut suppresses the cross section by above 10% in the 
whole rapidity range. For e=0.1 and R = 1 the suppres- 
sion is 17-23% at rapidities —1.5 < ?7 7 < 4. This large 
effect is not too sensitive to the value of e: changing the 
value by a factor of 2 from e = 0.1 to e = 0.2 or to 
e = 0.05 varies the results for isolated photon by about 
4% (fig. 0a). The dependence on R is stronger but also 
not very large: when changing R value by a factor of 2 
from 1 to 0.5 the results increase by about 7% (fig. 0b). 

The suppression due to the isolation imposed on the 
photon is presented in Table 1 (the second line) for indi- 
vidual contributions and for the total cross section. As 
expected, the cross section for fragmentation processes 
(i.e. with resolved final photon) is strongly suppressed: 
after isolation it is lowered by a factor of 5. At the same 
time the QCD corrections to the Born diagram increase 
significantly, i.e. the contribution to the subtraction cross 
section, dcr su b, due to this corrections is negative. The 
isolation restrictions do not modify contributions of other 
subprocesses since they involve photons isolated by def- 
inition. The subtraction cross section, being a sum of 
negative QCD corrections and fragmentation contribu- 
tions, is of course positive and the total cross section for 
isolated final photon is lower, by 20%, than for the non- 
isolated one. 

In the following we keep R=l and e=0.1, standard 
values used in both theoretical and experimental analy- 
ses. 

5.2 Other experimental cuts 

In order to compare the results with data we consider 
other cuts imposed by the ZEIIS group on prompt photon 
events at the HERA colliderE3. The influence of the lim- 
ited energy range, 0.2 < y < 0.9, is shown in fig. [| The 
cross section is strongly reduced, by 30-85%, in the pos- 
itive rapidity region. At negative rapidities the change 
due to the y-cut is weaker: 5-10% at —1.2 < ry 7 < —0.4 
and 10-30% at other negative rapidities. Separately we 
show the results obtained without the box subprocess 
(fig. |8|). The box diagram contributes mainly in the ra- 
pidity region between -1 and 3. After imposing the y-cut 
its contribution is important in narrower region from -1 
to 1. The influence of the y-cut can be read also from 
Table 1 (the third line). One sees e.g. that the Born 
contribution is reduced 3.5 times, while other ones are 
suppressed less, roughly by a factor of 2. 

The results obtained for the isolated photon with the 
y-cut and in addition with the cut on the final photon 
rapidity, —0.7 < rf < 0.9, are presented in the last line 
of Table 1. The restriction on rf decreases the contri- 
butions of all subprocesses approximately by a factor of 
2 (except for the double resolved contribution reduced 
almost 3 times). 
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The role of various experimental cuts is illustrated 
also in fig. |^, this time for the x 1 distribution. In partic- 
ular we see that the isolation and the energy cut reduce 
considerably the contributions from large and medium 
x 1 , while the contributions from x-y below 0.1 are reduced 
less. On the other hand, the small x 1 contributions are 
strongly, by two orders of magnitude, diminished by the 
photon rapidity cut. This shows that measurements at 
the central rf region (—0.7 < if < 0.9) are not too sen- 
sitive to the small x 7 values in the photon. 

When calculating the QCD corrections to the Born 
process in the subtraction term da su b we used the small-<5 
approximation described in sec. |^. Because these correc- 
tions give less than 10% of the cross section for the iso- 
lated photon production with various cuts (see the third 
column in Table 1), we expect that the error resulting 
from using the approximations is small, though we use 
in fact not a small value of 5 (5 — R/coshrj^, R = 1). ^| 

5.3 The comparison with data 

Two types of the final state were measured in the ZEUS 
experiment: 1) an isolated photon with —0.7 <rf< 0.9 
and 5 < pr < 10 GeV; 2) an isolated photon plus jet 
with the photon rapidity and transverse momentum as 
above, the jet rapidity in the range —1.5 < 77 jet < 1.8, 
and the jet transverse momentum > 5 GeV. 

We compare our NLO predictions with the ZEUS 
data from the first type of measurements!! 2 !. In fig. |lO| a 
the comparison is made for the transverse momentum 
distribution for various Nf. Although the predictions 
tend to lie slightly below the data a satisfactory agree- 
ment is obtained for Nf = 4. Note a large difference 
between the results for Nf—A and 3 due to the fourth 
power of electric charge characterizing processes with two 
photons. We observe a very small contribution from the 
bottom quark (for Nf=5). The predictions are obtained 
in massless quark scheme and may overestimate the pro- 
duction rate. 

Similar comparison of the NLO results with the data, 
now for the rapidity distribution, is shown in fig. [l0|b. 
A good description of the data is obtained for Nf=4 
and Nf—5 in the rapidity region 0.1 < ?y 7 < 0.9. For 
—0.7 < r/j < 0.1 our predictions lie mostly below the ex- 
perimental points. This disagreement between predicted 
and measured cross sections is observed also for other 
theoretical calculations (LG) and for Monte Carlo sim- 

c Calculations fp« the prompt photon production in ep and 
hadron-hadron E3 collisions were performed using space slicing 
method without the small 8 assumption. Comparison of such re- 
sults with predictions obtained in (discussed here) approximated 
way showed that the small 8 approximation is an accurate analytic 
technique fpfl including isolation effects in NLO calculations (also 
for R = 1) 15. 



ulations In fig. [ToJd we present separately an effect 
due to the box subprocess (for Nf =4). It is clear that 
the box term enhances considerably the cross section in 
the measured rapidity region. Its contribution to the 
integrated cross section is equal to 9.6%. The double 
resolved photon contribution is also sizeable, although 
roughly two times smaller than the box one, see Table 1 
(fourth line). Both these [a 2 ] contributions improve de- 
scription of the data. 

The predictions obtained using threeJifferent_NLO 
parton densities in the photon (ACFGP 0, GRVtS and 
GSl3) are presented for Nf = 4 in fig. |TI| a (Q = pr) and 
in fig. [ll]b (Q = 2p T ) together with the ZEUS datalii 
The results based on ACFGP and GRV parametrizations 
differ by less than 4% at rapidities ?y 7 < 1 (at higher t/ 7 
the difference is bigger), and both give good description 
of the data in the rapidity range 0.1 < r/y < 0.9 (for 
Q = Pt and Q — 2p T ). For —0.7 < ?y 7 < 0.1 none of 
the predictions is in agreement with the measured cross 
section. 

For Q=pr (fig- |!]a) the GS distribution leads to 
results considerably below ones obtained using ACFGP 
and GRV densities, especially in the rapidity region from 
roughly -1 to 1. This difference between the GS and other 
considered herein parton parametrizations is mainly due 
to their different treatment of the charm quark in the 
photon. In the GS approach the charm quark is absent 
for Q 2 below 50 GeV 2 . Since we take 5 < Q = p T < 10 
GeV, and the most important contribution to the cross 
section arises from the lower pt region (see fig. |To|a) , the 
Q 2 value usually lies below the GS charm quark thresh- 
old. As a consequence, predictions based on GS have 
strongly suppressed the contribution of subprocesses in- 
volving charm from the photon - contrary to GRV and 
ACFGP predictions where the charm threshold is at 
lower Q 2 . 

The above explanation of differences between cross 
sections involving GS and both GRV and ACFGP par- 
ton densities is insufficient for higher rapidities, ?7 7 > 2. 
Here the differences between the results based on partic- 
ular photon parametrizations are bigger, especially when 
comparing predictions obtained using GRV and ACFGP 
ones (not shown). This is due to large differences be- 
tween used parton densities at low £ 7 , which is probed 
at the high rapidity region. 

All the considered parton distributions give similar 
description of the data when the scale is changed to 

d The ZEUS Collaboration has presented recently an analysis Q 
in which an intrinsic transverse momentum of partons in the pro- 
ton, kx, was introduced in the PYTHIA 6.f generator in order 
to improve agreement between the data and Monte Carlo predic- 
tions for an isolated photon plus jet photoproduction. The data, 
selected with rr 7 > 0.9, are consistent with the predictions for 
< k T >= 1 .69 ± 0.181°;^ GeV. 
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Q = 2pT : see fig. p"T|b. Here the calculation corresponds 
to Q 2 which is always above 50 GeV 2 and the charm den- 
sity in the GS parametrization is non-zero, as in other 
parametrizations . 

In fig. [12] our predictions are compared to the ZEUS 
data divided into three ranges of y. This allows to estab- 
lish that the above discussed discrepancy between the 
data and the predictions for rp < 0.1 is coming mainly 
from the low y region, 0.2 < y < 0.32. In the high y 
region, 0.5 < y < 0.9, a good agreement is obtained. 

We have also studied the dependence of our re- 
sults on the choice of the parton distributions in the 
proton and parton fragmentation into the photon (not 
shown). Cioss sections calculated using GRV El, MRST 
(set ft08a)E9 and CTEQ4M0NLO parton parametriza- 
tions for the proton vary among one another by 4 to 7% at 
negative rapidities and less than 4% at positive rapidity 
values. Results for the isolated final photon are also not 
too sensitive to the fragmentation function. For rapid- 
ity ranging from -1 to 4 the cross section obtained with 
DO LOa fragmentation function is 2-=- 3.5% lower than 
the cross section based on GRV NLOtiil parametrization. 
Only at minimal (ry 7 < —1) and maximal (4 < 7y 7 ) ra- 
pidity values this difference is larger, being at a level of 
3.5 - 8%. 

6 The theoretical uncertainties of the results 
and comparison with other NLO predictions 

As we already mentioned the predictions are obtained in 
massless quark scheme and may overestimate the produc- 
tion rate. An improved treatment of the charm quark, 
especially in the box contribution which is particularly 
sensitive to the change from Nf — 3 to Nf — 4, would 
be needed. However we do not expect that this improve- 
ment would change qualitatively our results. 

We now discuss the theoretical uncertainties of our 
predictions related to the perturbative expansion. 

6.1 The dependence on the Q scale 

In order to estimate the contribution due to missing 
higher order terms, the influence of the choice of the 
Q scale is studied for the ry 7 distribution. In fig. |l3| a 
the results obtained using GRV densities with and with- 
out the y-cut are shown. When changing Q from px to 
2pT (pr/2) the cross section increases (decreases) at ra- 
pidities below ~ 1 and decreases (increases) at higher 
rapidity values. Only at high rapidities (where the cross 
section is small), 7y 7 > 3, the dependence on the choice of 
the scale is strong, above 10% (up to 20-30% at 7y 7 « 5). 
In the wide kinematical region, —2 < rj 7 < 2, the rel- 
ative differences between results (with and without the 
y-cut) for Q ~ pt and results for Q = 2px or pr/2 are 



small and do not exceed 6%. Around the maximum of 
the cross section at rapidities — 1 < ry 7 < these dif- 
ferences are 4-6%. This small sensitivity of the results 
to the change of the scale is important since it indicates 
that the contribution from neglected NNLO and higher 
order terms is not significant. 

Note that individual contributions are strongly de- 
pendent on the choice of Q, e.g. results for the single 
resolved processes vary by ±10-20% at rapidities r] 7 < 1. 
Results are much more stable only when the sum of re- 
solved processes and QCD corrections is considered. 

In fig. [l3] b we present NLO results for various Q with 
and without the y-cut, however this time with no box 
contribution. At rapidities rj 7 < 1 the uncertainty due 
to the choice of the renormalization scale is about two 
times higher than for the cross section with included box 
diagram, so the box contribution (~ [a^]) seems to stabi- 
lize the NLO prediction. At rapidities r) 7 > 2 the relative 
dependence on the choice of the scale is similar for the 
cross section with and without the box term. 

6.2 The comparison of NLO and LL predictions 

In the present calculation we include in LL accuracy the 
single and double resolved photon processes as well as 
the bjoxdiagram in addition to the Born contribution, see 
alsotl cm (although this is not fully consistent with the 
discussion in sec. |J). The cross section for the jp — > 7 A 
scattering in LL accuracy is obtained by convolution of 
partonic cross sections with relevant LO parton densities. 

In fig. [[4] we show the LL prediction for the isolated 
7 photoproduction (dotted line) together with NLO pre- 
dictions (solid line) and the Born contribution only (dot- 
dashed line). The highest differences between the LL 
and NLO cross sections are seen at the rapidity range 
-0.5 < 7y 7 < 2.5 where the LL results lie 10-20% below 
the NLO ones (fig. |i~4|a). For the pr distribution this 
difference is 10-14% in the whole presented range of the 
transverse momentum, 4 < px < 20 GeV (fig. pT|b). 

We think that the observed difference between NLO 
and LL results together with the weak dependence on 
the Q scale discussed in sec. |6.l| indicates reliability of 
the calculation. 

6.3 The comparison with other NLO results 

As we discussed in Sec. 0, our NLO calculation of the 
DIC process diffkrsJxom the "l/a s "-type NLO analysis 
presented in ref.lij l3, by set of diagrams included in the 
calculation. We do not take into account a s corrections 
to the single and double resolved processes, which are be- 
yond the NLO accuracy in our approach. On thepOther 
hand, we include the box diagram neglected in Einl3. 
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(The double resolved subprocesses are included in both 
analyses.) 

We compare our results and the results of the LG 
calculation t3 (using Nf=A and Q = px) for the isolated 
final photon (i?=l r £ | =0.1) in the kinematical range as in 
the ZEUS analysisti (i.e. for -0.7 < rp < 0.9 and 0.2 < 
y < 0.9). First we use the GRV photon parton densities. 
The LG predictions for da/dpx cross section are about 
20% higher than ours in the presented range of transverse 
momentum, 4 < px < 20 GeV. For da/drf cross section 
(with 5 < px < 10 GeV) the biggest differences are at 
r/ 7 =0.9 where the LG results are about 35% higher. The 
differences decrease towards negative rapidity values and 
are negligible at —0.7 < rf < —0.5. For y range limited 
to low values only, 0.2 < y < 0.32, the LG cross section 
is higher than ours by up to 20% at positive f? 7 , while at 
negative r; 7 it is lower by up to 10%. For large y values, 
0.5 < y < 0.9, where our predictions agree with data, the 
LG results are higher than ours by up to 80% (at r\ n = 
0.9). 

As it was already discussed, for Q = p T the GS pho- 
ton distributions lead to results lower by 11-14% than 
those obtained with GRV densities at rapidities betase 
-1 and 1 (see sec. In calculations presented inEiHi 
this difference is even twice larger. 

The LG predictions obtained using the GS 
parametrization lie up to 20% below ours (also based on 
GS distributions, with Q = px and Nf = 4) at rapidities 
—0.7 < 7/ 7 < 0.2, and th^y_are higher than ours by up to 
30% for 0.2 < ?7 7 < 0.9 



7 The LL prediction for 7+ jet photoproduction 

The ZEUS Collaboration has also analyzed the prompt 
photon photopmdjicjLion in which in addition a hadron 
jet is measuredlij'OEl In fig. |lj we show the LL predic- 
tion for the isolated 7+ jet final state ^ together with the 
predictions for the 7 alone. The following jet rapidity and 
transverse momentum are assumed: —1.5 < rjjet < 1-8 
and p!^ 1 > 5 GeV, respectively. This additional cuts for 
the final state imposed on jets decrease the cross section, 
especially at high rapidities. The LL predictions for the 
7+jei are lower than those for the 7 production by 5-10% 
at negative rapidities and by 10-80% at positive rapidity 
values. The difference between both LL results is about 
10% in wide range of transverse momenta, ~ 6 < px < 20 
GeV, and only for the lower px region, 4 < px < 6 GeV, 
it is higher (13-23%). 



e The NLO calculation for 7+j'et photoproduction will be discussed 
in our next paper. 



8 Summary 

Results of the NLO calculation, with NNLO contribu- 
tions from double resolved photon processes and box di- 
agram, for the isolated 7 production in the DIC process 
at HERA are presented Q The role oLthe kinematical 
cuts used in the ZEUS measurement E£l are studied in 
detail. 

The results obtained using GRV parametrizations 
agree with the data in shape and normalization for px 
distribution. For rf distribution a good description of 
the data is obtained for rf > 0.1, while for if < 0.1 the 
data usually lie above the predictions. This discrepancy 
arises mainly from the low y region, 0.2 < y < 0.32. 
The beyond NLO terms, especially a box contribution, 
improve the description of the data. 

We have studied the theoretical uncertainty of results 
due to the choice of the renormalization/factorization 
scale: Q — px/2,Pt,2pt- At high rapidities rj 7 > 3, 
where the cross section is small, the uncertainty is 10- 
30%. In a wide range of rapidities, —2 < ?] y < 2, the 
dependence on the Q scale is small, below 6%. Since we 
include some NNLO diagrams in our NLO calculation, 
this stability of the predictions versus the change of the 
scale is especially important. The week dependence on 
the Q scale, and not large differences between LL and 
NLO predictions (below 20%) allows to conclude that 
theoretical uncertainties of our NLO calculations for an 
isolated photon production in the DIC process at HERA 
are relatively small. 

We compared our results with the LG ones based on 
different set of subprocesses. The cross section da/dpx 
obtained by LG is about 20% higher than ours (for GRV 
photonic parton distributions). For the cross section 
da/di] 7 this difference is up to 35% at ?7 7 = 0.9. The 
highest differences are present for high y values only, 
0.5 < y < 0.9, where on the other hand our predic- 
tions are in agreement with the data. At low y range, 
0.2 < y < 0.32, differences between both calculations 
are smaller and none of them describe the data well for 
rapidities below 0.1. 
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Appendix 

Here we present formulae for the subtraction term in the cross section for the production of the isolated photon (see 
sec. |^) in the jp — > -fX scattering. The corresponding expressions for the ep reaction can be obtained using the 
equivalent photon approximation (sec. ^). The subtraction term is the cross section for subprocesses in which the energy 
of hadrons inside the cone of radius R around the final photon is higher than the energy of the photon multiplied by a 
small parameter e, 

J2^h>eE y . (12) 

h 

The cone is defined in the rapidity and azimuthal angle plane: 

y/(Vh-rh) 2 + (<l>h-^) 2 <R. (13) 

Below we use the variables v and w defined in the following way: 

t —it 

v = 1 + - , w = K, 

s s + t 

where s, t and u are the Mandelstam variables for the partonic subprocess, 

s — ijx^xSq'p ■ t — yx^ zTq^p . u — x z\J r g>p ^ 
and S e p, T ep and U ep stand for the Mandelstam variables for ep — > ejX reaction, 

Sep — ^PePp , T e p — 2p e p 7 , Uep — ^PpP^ • 

The p e , p p and p 1 denote the initial electron and proton four- momenta and the four-momentum of the final photon, 
respectively. The fractional momenta y, x 1 , x and z are defined in sees. 



The subtraction cross section consists from two contributions which arise from subprocesses involving fragmentation 
function and from a s corrections to the Born process: 



2.2 and 



3 ip-nx d 3 cr1 P ~ i ' lX d 3 a^^ x 



7 d 3 Pl 



_ £. f rag ^ ^ 

7 <i 3 p 7 7 d 3 p 1 



J^J ™ ~ Sub J- J- " ' n,r 



(14) 



where 



frag 



d 3 p. 



= E E / — / dxf b/ p{x,Q 2 )D l/c {z,Q 2 )E, 



b=g,q,qc=g,q,q 



d^cjfb^cd 

d 3 p- ) 



E E E / — dx J dxf ah {x 1 ,Q 2 )f b /p{x,Q 2 )D l/c {z,Q 2 )E r 
= t = „ „Jo z Jo Jo 



d 3 a a 



b — >cd 



a=g,q,qb=g,q,qc=g,q,q 



d 3 p~ 



(15) 



and 



J3 7P~>7* ^ /"I 



d 3 P-) 



v(l — w) 



j3_7<?i-^7<?i+S j3 7<?i^79+9; 

Q 2 )£ 7 s g + 4 /p (x, Q 2 )£ 7 ^ 



d 3 p 7 



+ / ff/p (z,<2 2 )i? 7 ^^ 



3^79^79i+9i 



d 3 p 7 



(16) 
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with 



x 



-yT ep 
yS ep -\-U e 



The contribution (|15|) comes from 2 — > 2 single and double resolved subprocessesJrL which the final c-parton decays into 
the photon. Here the calculations are standard, as for non-isolated photon caseclnj. The condition ( |l2| ) is included via 
the upper limit of the integration over z: z = Ey/E c < 1/(1 + e) (the hadron remnant of the c-parton takes the energy 
(l/z — \)E 1 > eEj). The lower limit of the integration over the fractional momenta z, x 7 and 



Eh — E r 



E 1 , so E h 



x is formally zero but in fact this zero-value is inaccessible due to the delta function, 8{yx 1 xS ep + yx 1 zT ep + xzU ep ), 
the partonic cross sections for 2 — ► 2 subprocesses. 



The contribution ( |jq ) describes the a s corrections to the Born process. The diagrams involving the virtual gluon ex- 
change do not contribute to the subtraction term, and only 2 — > 3 processes are included. In these processes the photon 
and two partons are produced. One of the partons enters the cone of radius R around the photon, and its energy should 
be higher than the photon's energy multiplied by e. To fulfill this condition the integration is performed over x values 
for which v(l — w)/(l — v + vw) > e. There are three types of such processes: 

• 7(7 — > 7(7 + g (with a quark inside the cone) , 

• 19 ~^ 19 + 1 (with a gluon inside the cone), 

• 19 ~~ * 11 + ? (with a quark inside the cone) . 



Below we present our analytical results for the a s corrections contributing to the subtraction term. The results are 
obtained with the assumption that the angle between the photon and the parton inside the cone is small. The quark 
masses are neglected. All collinear singularities are shifted into the fragmentation functions D (infrared singularities 
do not appear in this calculations). The px stands for the transverse momentum of the final photon. The partonic cross 
sections in (hq) are given by following expressions: 



E~ 



d A <7 



sub 



d 3 P~, 



A a 2 



ITS 



(1 - v + vw) 2 + (1 - vf - 2 + (R PT f 
(1 — v + vw)(l — v) * s 



(vwY 



(1 



vw)(l 



(17) 



J3_79i^7ff+9i 
P " "sub 

7 d3p 7 



^ a em a s 4 
3 7TS 2 qi 



(Rpr) (1 - v) (1 - v + vw) 2 + (vw) 2 ri . , 4 4 ,, , 4l 

. ±jz — — ^ 7- ' — . / 1 i + (i - v + vwf + v*(i - wr] , (18) 

s (1 - v + vw) 5 v(l ~ w){vw) 2 L v Jl 



E. 



<Pp^ 



1 o.„ m a 



(vw) 2 + (1 - vf 
vw(l — v) 



P(Q 2 



{Rpr? (1 



vw) 



(vw) 2 (l — v) 2 



(19) 



where 



p{Q2) = l + v 2 {l-w) 2 ^ (R 2 p 2 T v 2 {l - w ^ 2 
(1 — v + vw) 2 



(20) 
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Table 1: The cross section for non-isolated and isolated final photon, isolated photon with 0.2 < y < 0.9, and isolated photon with 
0-2 < y < 0.9 and -0.7 < rp < 0.9. 



[Pb] 


total 


Born 


0(a s ) 


box 


single resolved 


single resolved 


double resolved 












initial 7 


final 7 




non-isolated 


226.2 


82.1 


8.7 


13.9 


54.7 


24.6 


42.2 






(36.3%) 


(3.8%) 


(6.1%) 


(24.2%) 


(10.9%) 


(18.7%) 


isolated 


180.4 


82.1 


15.2 


13.9 


54.7 


5.12 


9.37 






(45.5%) 


(8.4%) 


(7.7%) 


(30.3%) 


(2.8%) 


(5.2%) 


isolated 


72.33 


23.6 


6.33 


6.54 


28.2 


2.34 


5.29 


y cut 




(32.6%) 


(8.8%) 


(9.0%) 


(39.0%) 


(3.2%) 


(7.3%) 


isolated 


35.36 


13.6 


3.32 


3.41 


11.9 


1.21 


1.92 


y, rj 7 cuts 




(38.5%) 


(9.4%) 


(9.6%) 


(33.7%) 


(3.4%) 


(5.4%) 



ep — >■ e^fX 

1000 1 1 1 1 1 1 1 — 1 — 1 




Figure 6: The final photon py-dcpendence of the cross section dcr/dpT for non-isolated 7 photoproduction (solid line). The Born contribution 
is shown separately (dash-dotted line). 
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ep — > e^X 



ep 



e^X 




b 



rr — i 1 1 — 

— ^ \ non-isolated 7 




Figure 7: The differential cross section da/drj^ as a function of the photon rapidity r\ 1 . a) The results for non- isolated photon (dashed line) 
and isolated photon with R = 1 and e = 0.05 (dotted line), 0.1 (solid line) and 0.2 (dot-dashed line), b) The results for non-isolated photon 
(dashed line) and isolated photon with e = 0.1 and R = 0.1 (dot-dashed line), 0.5 (dotted line) and 1 (solid line). 
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1 1 
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- 1 ^ \y 
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/]/ without box 




\ \ \ \ 
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Figure 8: The difTercntial cross section da/drjj for isolated 7 (e = 0.1, R = 1) as a function of the photon rapidity ?? 7 with (solid lines) and 
without (dashed lines) the box contribution. The results arc obtained with imposed y cut (0.2 < y < 0.9) and without this cut. 
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100 |r Q = pt 

N S =A 
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0.01 
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_L 
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Figure 9: The cross section in m 7 bins of the length 0.1. The results for non-isolated 7 integrated over the whole range of y and r/j are shown 
with the dashed line. The solid line represents results integrated over the whole range of y and for isolated 7 with e = 0.1 and R = 1. 
Results with additional cuts in the isolated 7 cross section are shown with: dotted line (0.2 <y< 0.9) and dot-dashed line (0.2 < y < 0.9, 
-0.7 < »? 7 < 0.9). 



ep — > cyX ep — > ejX 




4 6 8 10 12 14 16 18 20 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

p T [GeV] % 



Figure 10: The results for isolated]— « with various numbers of active massless flavors: Nf = 3 (dashed lines), 4 (solid lines) and 5 (dotted 
lines), compared to the ZEUS dataE3. a) The differential cross section da/dpx as a function of the photon transverse momentum, b) The 
differential cross section da/drj^ as a function of the photon rapidity the result without the box contribution is also shown for Nf = 4 
(dot-dashed line). 
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ep -> e^X ep — s> ejX 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

r? 7 77 7 



Figure 11: The differential cross section da/dr/-/ for isolated 7 pp a function of the photon rapidity ?j 7 compared to the ZEUS datat3. Three 
different NLO photon parton distributions are used: ACFGPH (dotted linrA^GRVEl (solid line) and G_S EZI (dashed line). The GRV NLO 
parton distributions in the proton£3 and parton fragmentation into photon £3 are used, a) Q = px- b) Q = 2px- 
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Figure 12: The results for three ranges of y: 0.2 < y < 0.32, 0.32 < y < 0.5 and 0.5 < y < 0.9, compared to the ZEUS data 



Ei. 



15 



ep — > e^fX ep — ► e^X 




Figure 13: The differential cross section dcr/d^ for isolated 7 photoproduction as a function of the photon rapidity -q-y with (a) and without 
(b) the box contribution. Three different values of the Q scale are assumed: Q = (dashed lines), Q = Pt (solid lines) and Q = 2px 

(dotted lines). The results are obtained with imposed y cut (0.2 < y < 0.9) and without this cut. 




Figure 14: The differential cross sections dcr/dr)^ (a) and da/dpT (b): NLO (solid line) and LL (dotted line) results for isolated 7, and LL 
predictions for isolated 7 + jet photoproduction (dashed line); the jet rapidity is assumed in the range: —1.5 < ry e t < 1.8. The dot-dashed 
lines show the Born contribution to the NLO cross section for isolated 7. GRV NLO (LO) parton densities in the proton and photon were 
applied in NLO (LL) calculations. 
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